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■^ We extend our earlier model of the small scale structure of cosmic string networks 

^} through an improved treatment of the separation of long and short scales. We find that 

^^ the production of small loops (at the gravitational radiation scale) is a robust feature 

] of string networks, in addition to a population of loops near the horizon scale. We 

T— I obtain quantitative agreement with the scaling of loop production functions as found 

^^ in simulations by two groups. 

> 

X 



^dubathOkitp . ucsb . edu 
^ j oepOkitp . ucsb . edu 
■^jrochaSphysics .ucsb. edu 



Estimates of the size of cosmic string loops are uncertain to many orders of magnitudes, 
ranging from near-horizon scales [I1I21E1I1] down to the gravitational radiation length scale ^ 
El [71 [8], and even to microscopic scales [9]. Even when the properties of the string are 
specified, so that the problem is completely well-posed, this uncertainty arises from the 
inability to solve the equations for the evolution of the network. Numerical simulations have 
been unable to reach the necessary dynamical ranges of length and time, while analytic 
methods have difficulty with the nonlinearities of the problem. 

There has recently been progress in both approaches. By using a doubling trick, recent 
simulations have been able to reach a much larger dynamical range [21 H] • By focusing on the 
properties of the network well below the horizon scale, it has been possible to formulate an 
analytic model that captures the quantitative properties of the short distance structure [10] . 
In this paper, we will argue that these two recent advances point to a common picture. 

The new simulations, taken at face value, show two peaks in the loop production func- 
tion [31 H] . One is at a surprisingly large scale, around a tenth of the horizon length, while 
the other is near the UV cutoff. The large loops have strong gravitational wave signatures, 
and so lead to stronger bounds and enhanced prospects for future sensitivity [HI [121 [T3] - 
However, the total length going into loops is fixed by energy conservation, so any string 
that goes into small loops reduces this effect [H]. The simulations at their current dynamic 
ranges show roughly 90% of the string going into small loops, but the authors of that work 
conjecture that these loops are transient, so that in final attractor solution essentially all 
string will go into long loops. 

The analytic model begins by considering the small scale structure on long strings [TO]. 
It is found that the string approaches fractal dimension 1 at short distances, with a deviation 
that goes as a power of the separation /, 



df = l + 0{l/t) 



2x 



The exponent x is determined in terms of the rate of expansion of the universe, a{t) oc f^ 
and the mean velocity squared f ^ in the network, taken from simulations: 

^ l-z/(l-2i;2) ^ ' 

The relevant values are 

radiation era : u = - , v'^ = 0.41 , x = 0.10 , 

matter era : z/ = - , t>^ = 0.35 , x = 0-25 . (3) 

o 

The predicted two-point function appears to agree with simulations [21 [15]. Taking this 
structure as input, we can then calculate the rate at which self-intersections of the long 
string take place to produce a loop, as a function of loop size. One finds that, even though 
the string is becoming straight at short distance, the loop production function (that is, the 
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rate of loop production weighted by loop length) diverges at small lengths [TO]. This is 
because the exponent x is rather small in both eras, meaning that the approach to fractal 
dimension 1 is slow: the critical value for small loop production is x = O.Sq 

The divergence in the loop production function must be cut off, because the total rate 
at which string length goes into loops is fixed by energy conservation. In the approach of 
Ref. [To], this condition saturated at a length scale not too far below the horizon. It was then 
argued that smaller loops would be produced by extensive fragmentation, but this appeared 
to be highly nonlinear and not amenable to analytic calculation. In the present paper we 
make an essential improvement on the previous model, which leads to a very different, and 
simpler, picture. 

On the time scale of small loop formation, we can ignore the expansion of the universe. 
The string configuration is described by the left- and right-moving unit tangent vectors 

(p+(M),p_(f)) = 2(9„x(M,t;),9^,x(M,f)) , {u,v) = {t + a,t- a) . (4) 

The calculation of the loop production rate requires separating the p-t into a long-distance 
piece and a short-distance piece. In Ref. [10] the long-distance part was simply averaged over 
the unit sphere, with a weight factor chosen to give the correct mean value of p+ ■ p_. In the 
present work we will take a long-distance part with some given classical (m, f )-dependence, 
and then average over the short distance part with a weighting corresponding to Eq. ([I]). In 
the end we should average over an ensemble of long-distance configurations, but the result 
will be largely independent of this. The inclusion of the time-dependence will reveal that 
loop formation leads to a strong correction to the distribution of long-distance configurations, 
which was not accounted for in the earlier approach [TO]F^ 

A typical configuration is shown in Fig. 1. The long-distance configuration is in bold, 
and the total configuration is dashed. Note that, while the string configuration approaches 
fractal dimension 1 at short distance, the graph of its tangent has large fractal dimension 

In the figure, we have shown a long distance configuration in which the curves p+(m) and 
p„(f) cross, meaning that there is a cusp in the spacetime evolution of the string [T^. It 
is evident that the effect of the short distance structure is to turn this one large cusp into 
many small cusps. The formation of small loops is closely associated with the cusps. The 
condition for a cusp is crossing of the two curves at a point p+('u) = P-('y). The condition 
for self-intersection of the string is 

L+(n,/) = L_(t;,/) , 

L+{u,l)= du'p+{u'), L_{v,l)= dv'p4v'); (5) 

Ju-l/2 Jv-l/2 



^In Ref. [2 the loop production function shows no evidence of this divergence but this is expected: in 
those simulations the short distance structure seems to change below a certain length scale in such a way 
that X exceeds the critical value. This could be an artifact of the finiteness of their dynamic range. 

^These same issues were faced in the pioneering work 16J , in a somewhat different framework. 
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Figure 1: The functions p± written as a fixed long-distance piece (bold) plus a random 
short-distance part. The values of u and v increase along the arrows, with the classical cusp 
defined to lie at {u,v) = (0,0). 



when this is satisfied the segment of length /, centered at {u,v), breaks off as a loop. For 
small / the functions p-t will have small variation along the segment, so small loops will form 
where these are approximately equal [IHl [IS] . 

Following this discussion, it is clear that at early times, when u and/or v are sufficiently 
negative, the curves will be outside the cusp region and the production of small loops will 
be negligible. As the cusp begins to form, production of small loops begins to increase. We 
will find that loops of all scales begin to form at the same time, extending all the way down 
to the gravitational radiation cutoff. This is in marked contrast to the cascade from large 
to small loops, anticipated in Ref. ^U\. Thus we are led to conclude that the production of 
loops down to the gravitational radiation scale is inevitable, and not a transient effect. 

The production of large loops is outside our model, but these also seem inevitable. In 
the course of the evolution of the network it must happen at least occasionally that a large 
loop will form and survive reconnection. We must leave it to the simulations to determine 
the relative weight and average size of these loops, but we can understand why they do not 
fragment entirely into much smaller loops. From Fig. 1 we expect that if there are any cusps 
on these loops then the regions containing them will dissolve into small loops. However, if 
there are parts of the loop where p+ is sufficiently far from the p_ curve, and vice versa, 
then these will survive. Thus we expect that if the functions p± are plotted, there will be 
significant gaps in each, corresponding to a large kink, and the surviving curves will avoid 
each other. It seems plausible that the total amount of string surviving in the long loops is 
a small fraction of the total, as seems to be true with the current dynamical range [31 0], but 



it is important to study this further. 

The loop production cuts off when the probabihty for a given bit of string to be incor- 
porated in a loop reaches 1. We will see that this occurs long before the classical cusp point 
{u, v) = (0, 0) is reached. Thus, whereas in Ref. |10] the region where the long-distance parts 
of p+{u) and P-{v) were parallel gave the dominant contribution, but in a correct treatment 
it is simply absent. This explains why this earlier treatment gave a normalization for the 
loop production function that was much too large, though it did seem to find the correct 
functional form. 

In the remainder of this paper we will carry out the calculation just described. As 
further support for our picture, we will obtain a prediction for the small loop production 
function that fits well with the simulations by two groups. The fit actually works better than 
anticipated, as there is no sign that it is distorted by further fragmentation; we speculate as 
to why this might be. 

To carry out the calculation of the rate, we separate the left- and right-moving unit 
tangent vectors into their long- and short- distance parts, 

P± = 1± + s± . (6) 

Since the production of small loops takes place near the cusp, we take the simplest cusp 
form for the long-distance structure, 

1+ = z(i - vyy/' + v+n ^ z + v+u - ^vy , 

1_ = z(l - 1/_V)^/^ + V^v ^ z + V^v - -K_V (7) 

with V+ ■ z = V_ ■ z = 0. We are expanding to order u"^, v"^ . The slopes V-t have units of 
inverse length, and are reciprocal to the size of the cusp. They will be of order the inverse 
correlation length, somewhat less than the horizon length. 

Adding in the short distance structure gives, to the same order, 

p± = z + y± - 2^± ' 

y+(M) = V+m + w+(m) , y_(f) = V_f + w_(t;) , (8) 



with 



{w,(.) .w,(.')) = -4c, £ |..*-"')j^^(e-l'l'<=w _ 



e 



'^ Re[{i\u-u'\+Lf>'-{i\u-u'\+lG^^f>'] , (9) 



t^^ cos TIX 



and similarly for w_; here c^ = 2 sin(7rx)r(l + 2%) and A is simply related to the normaliza- 
tion of the 0{l/tY^ correction in equation ([T]). It is approximately equal to 0.6 during both 



radiation and matter domination [10]. We have introduced a smooth short- distance cutoff 
at the scale Zqw ~ A{GfiY~^^H [S] to take account of the effects of gravitational radiation. 
The cutoff L represents the separation between the long-distance and short-distance parts 
of the configuration, taken around the correlation length. In the range /qw ^ \u — u'\ <^ L 
expression (|9| reduces to 



(w+(n) ■ w+iu')) ^ ^(L'2>^ -\u- uf ^)/t2x 



(10) 



with L'^^ = L^^/cosTTx; the u,u' dependence is as in Ref. [TO]. 

For reconnection probability unity, which we assume, the rate of loop formation is given 
by the rate of self-intersection. 



dJ\f = ((5^(L+(m,/) - L^{v,l)) |det J(m,i;,/)|) dudvdl 
where 'L± were defined above and the Jacobian is 



'11^ 



J(ii,f , /) 



P+(iz + //2)-p+(m-//2) 
p_(i; + //2)-p_(i;-//2) 
L |[p+(n + 1/2) + p+(« - 1/2) -v-{v + 1/2) -v-{v- 1/2)] 



(12) 



In our previous work [10] we were able to estimate the expectation value of the product (11 ) 
by the product of expectation values, but here there are strong correlations: when the delta- 
function is nonzero, the third row of the Jacobian is much smaller than its mean value. We 
therefore go to new variables, separating wj. into a piece constant on the segment and a 
piece with zero average on the segment. 



w+(m') = W+ + (jjj^{u) 
W 



w_ f 



UJ_[V 



W 



u+l/2 



du w^(u) 



u-l/2 
-^ (■v+l/2 



I 



dv'w_(v') . 



(13) 



v-l/2 



The variables uj± and W-t depend on the parameters u,v,l of the loop, but we leave this 
implicit. We will be interested in the loop production in the range /gw -C / <^ Iv, so we use 



the form (10). One then finds the two-point functions 



t2>^(a;+(u') ■ cc;+(?/)) = -A\u' - n'f ^ + /(«') + /{u") 
fHu;^iu')-W^) = -f{u'), 

t2x(w+ • W+) = ^L'2x + 0(/2x) , 



where 



fin) 



A 



{1/2 + u' - u)'^+i + {1/2 -u' + u) 



2x+l 



(2x + l)^ L 
is defined only for points on the segment, i.e. m— |<m'<m + |. 



•) 


(14) 




(15) 




(16) 


2X + 2J 


(17) 



Now let us express the rate in terms of these quantities. First, for the transverse parts 
of L-t we have 

L^ = /(W+ + V+m) , L^ = /(W_+V_f) . (18) 

The transverse part of the delta-function sets these equal, L^ = Lf = L-*-, and we use this 
in evaluating the remaining terms. For the z component. 



1 l'U+l/2 

Ll = l-- du' [V+n' + W+ + u+{v!)X 

^ Ju~l/2 



1 1 pu+l/2 

= I - -L^' - - / du' [oj4u') + V4u' - u)f , (19) 

^' ^ Ju~l/2 

where we have used the fact that the mean value of uj^ is zero. In the integrand, the 
fluctuations of u}+ are of order l^, while the classical term is of order /, so we can drop the 
latter for small loops. Thus, 



Ll-L^_ = -l r' ' du' [u;^{u')]' + I r' ' dv' [u,^{v')f = 0(1'^'^) . (20) 

^ Ju-l/2 -^ Jv-l/2 



''U+l/2 1 rv+l/2 

du [u:+{u')f + - / 

-1/2 ^ Jv-l/2 

For the transverse part of the first row of J, 



p^{u + 1/2) - p^{u - 1/2) = u}+{u + 1/2) - u+{u - 1/2) + V/ 

= u:+{u + l/2)-u+{u-l/2) + 0{l) . (21) 

Similarly in the second and third rows, in the transverse terms we can replace p± with a;±, 
after imposing L:[: = L:*:. In the z component of the first row, 

p\{u + 1/2) - pl(u - 1/2) = -yL^ ■ [^+{u + 1/2) - u:+{u - 1/2)] 

-\[[u^+{u + l/2)f - K(« - l/2)f) + 0(/i+^) . (22) 

The first term is of order l^ and the second of order Z^^, but the first actually drops out. 
In all three rows of J one finds the same pattern, J^ = — L"*" ■ J-*" + J^[^]. Thus the first 
term in the 2;-column is linearly dependent on the other two columns, and drops out in the 
determinant, leaving J^ with p± replaced by lo. 

We have now expressed the terms multiplying S'^(L^ — L^) all in terms of cj. The 



cross-correlation (15) is of order /^^, smaller than the geometric mean L'^l^ of the diagonal 



correlators (14, 16). We can therefore ignore it, giving 



{6'{L^ - L_) |det J|> = (<52(Li - L^)) x (<5(L; - LI) |det J|>p^_^ . (23) 

Thus we have achieved our aim of factorizing the expectation value. In the second term the 
delta-function scales as Z~^~^^, inversely to its argument, while the columns of J scale as 
1^,1^,1"^^, giving the same overall /"^"'"^^ scaling as in the previous work [10]. In fact, this 
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second term is identical to that in Ref. [TO]. It was evaluated there with the assumption that 
the short- distance structure is gaussian, giving a numerical coefficient of 2.2A in the radiation 
era and 1.3^ in the matter era. As we will discuss below, the gaussian approximation may 
not be valid here, but this does not alter the scaling with /. 
It remains to evaluate 



{S\Lt-L-)) = l 



CPq 
<fq 



giq-(L:|;-Li) 



^-^L'2xi2q2/4^2x_|_j;q.(v_|_u-V_»;) 

(27r)2^ 

^2xg-(V+w-V_t))2t2x/^L'2x 



(24) 



In all, 



dU = CL'-2^e-(^+"-^-'')'*'^/-^^"'^ dudv-^, (25) 

where (with the gaussian approximation) C = 0.7 in the radiation era and C = 0.4 in 
the matter era. These constants are actually independent of A so the normalization of the 
two-point function only enters the calculation of loop production through the exponential 
suppression factor. 

Now consider a left-moving point with given u. The total probability per unit v that this 
point be incorporated into a loop is 

— = C'L'-2Xe-(V+n-V_^,)2t2x/^L'2x f dl ^ 
dv J /2-2x 

Note that we have replaced J du —>■ I to count the loops containing the given point. This 
diverges at small / [10], because we have not yet taken into account the smoothing due to 
gravitational radiation in the two-point function (|9|. The smoothed two-point function is 



quadratic at m' — n ^ 0, corresponding to the form (10) at x = 1: at the shortest distances 



the divergence is gone (see also Ref. |20]). Thus we cut the integral off at Iqw to get 

— ^l^^Le . (27) 

At early times, where v is large, the probability of a loop containing u is small due to the 
gaussian. However, l^^^ ^ is large compared to the other dimensionful quantities, so as we 
integrate in v we soon reach total probability 1 for a range of values of u near the cusp: this 
portion of the string is removed by loop production (see Fig. [2]). Therefore, for each value 
of the M-coordinate, the maximum value of v that may still not be included in a loop is the 
solution (when it exists) of the equation Vivmax) = 1, where 
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Figure 2: The curve p+ is excised in the shaded region of the {u, f )-plane by loop formation. 
The cusp occurs at u = v = and the curve that dehneates the excised region is Vmax{u)- 
Both u and v are expressed in units of the cosmological time t. For this particular example 
the values Gfi = 10^^ and cos{9) = 2/7r were used, where 9 denotes the angle between V+ 
and V_. 



We see that the excised region depends on the the angle 9 between V_|_ and V_ for each 
cusp, not to mention the ubiquitous Gfi. To obtain the normalization of the loop production 
function would require the exact knowledge of the probability distribution for the angles 9 
as well as the density of cusps. One can certainly obtain normalizations in accordance with 
the simulations but there are too many free parameters so the answer is inconclusive. 

These are our main results. First, loops of all sizes down to the UV cutoff form simul- 
taneously, rather than in a cascade of fragmentation. Therefore the production of the small 



loops is a robust physical result. The integral of the probability (27) for a point on the 
string to break off becomes large long before the cusp f = is reached, so the overcounting 
in Ref. [10] is explained. However, this different treatment of the long-distance structure 



does not affect the short distance part of the rate (25), so the /-dependence is the same as 
before. 
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Figure 3: Mean values of x'^df{x) in the radiation era (upper bold curve) and matter era 
(lower bold curve, reduced by a factor of 10 for clarity), from Ref. [1]. The dashed lines show 
the slope predicted by our model. 



Fig. p^ shows the results of Ref. [3] for the average value of x^/(x), where x = l/t and 



3, ,2 



f{x) = 2t^7 



dl du dv 



(29) 



(here 7 is a dimensionless constant: the scaling value of the long string length per unit volume 
is 7^/t^). The simulations show a power law distribution above the UV cutoff, with exponent 
that match rather well with our model (dashed line) in both the matter and radiation eras. 
Ref. [7] also finds power law distributions for x'^f{x), with exponents 



Ref. [7] : a; ' (radiation 



our model : x ' (radiation) 



X 



-0.5 



-0.5 



(matter) , 
(matter) . 



(30) 



The simulations [7] and [1] thus agree quite well with each other, both matching our model 
well in the matter era and giving a slightly steeper distribution in the radiation era. 

We have found that fragmentation is not necessary to produce the small loops, but it 
still raises a puzzle. If we were to apply our same calculation to study fragmentation of 
these small loops, we would find again a strong tendency to fragment into loops at the UV 
cutoff scale, giving a sharp peak there rather than a power law distribution. We believe 
that the point is that the short distance structure on these small loops is not that typical 
of the long strings, but is dominated by two or more large kinks moving in each direction. 
One kink in each direction forms with the loop, and the remainder are present in the pre- 
existing distribution. We have noted in Ref. [10] that the tails of the small scale structure 
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distribution are not gaussian but are dominated by a single large kink (again, this does not 
alter the scaling with /). Thus these kinks are likely associated with the tail of the loop 
production function. If so, the functions p-t on the small loops are dominated by large 
jumps, from the kinks, and will be unlikely to have cusps where smaller loops can form. 

It may be interesting to consider a renormalization group (RG) analysis, where the sep- 
aration between the long-distance and short-distance configurations is taken at a running 
length scale L. Here we just note a few scaling properties. Taking the derivative of the two- 
point function, one finds that at scale L the rate of change of the tangent vectors is of order 
(L/t)^^^^. Thus the typical cusp velocities \± are increased by this factor, and the size of 
each cusp correspondingly decreased. At the same time, the lengths of the classical curves 
\± are increased by this factor, and so the number of intersections between them increases 
as {L / 1)"'^^'^^ . This gives another way to understand the loop production functions: it is 
plausible that the loop production per logarithmic scale IdM /dl scales as the number of cusps 
with L r^ I. This RG approach may be useful for understanding the spatial distribution of 
loop production, and correlations between loops. 

In conclusion, it appears that numerical and analytic methods are converging on a firm 
picture of the loop production, with one peak near the horizon scale and one near the 
UV cutoff, with perhaps the larger fraction of string in small loopsjj In Ref. [H] some 
observational consequences of this have been discussed. 

Note added: After the first version of this paper was posted, Ref. [23] appeared with 
a related title but different conclusions. Ref. |2S] contains several calculational errors that 
invalidate its conclusions. In particular, its quantity A3(A) is actually divergent, as follows 
immediately from the definition and the fact that 6 < 2 (note that the sign of c must also be 
corrected). This divergence has a simple origin. At any kink, a" will have a delta- function 
contribution, so the RMS value will necessarily diverge. Omitting this divergence essentially 
ignores the kinks, which of course gives an incorrect loop distribution. 
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case the cutoff may be at a smaller scale. 
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